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We consider the zero-temperature properties of the spin-half two-dimensional Shastry-Sutherland 
antiferromagnet by using a high-order coupled cluster method (CCM) treatment. We find that this 
model demonstrates various ground-state phases (Neel, magnetically disordered, orthogonal dimer), 
and we make predictions for the positions of the phase transition points. In particular, we find that 
orthogonal-dimer state becomes the ground state at J2/J1 ~ 1.477. For the critical point Jf/Ji 
where the semi-classical Neel order disappears we obtain a significantly lower value than J2/J1, 
namely, Jf/Ji in the range [1.14, 1.39]. We therefore conclude that an intermediate phase exists 
between the Neel and the dimer phases. An analysis of the energy of a competing spiral phase 
yields clear evidence that the spiral phase does not become the ground state for any value of J2. 
The intermediate phase is therefore magnetically disordered but may exhibit plaquette or columnar 
dimer ordering. 



I. INTRODUCTION 



The study of two-dimensional (2D) quantum mag- 
netism has attracted much experimental and theoreti- 
cal attention over many years. In 2D antiferromagnets 
at zero temperature the competition between interac- 
tions and quantum fluctuations is well balanced and one 
sees magnetic long-range order (LRO) as well as mag- 
netic disorder, dependent on details of the lattic o 1 ! 2 ^ 4 . 
In particular, frustration may lead to the breakdown of 
semi-classical Neel LRO in 2D quantum antiferromag- 
nets. Much research activity in this area has been focused 
on frustrated spin-half Heisenberg antiferromagnets on 
the square lattice, such as the J1-J2 model with com- 
peting antiferromagnetic nearest-neighbor J\ and next- 
nearest- neighbor J 2 bonds (see, e.g., Refsi 5 i 6 i 7 i 8 i 9 i 10 i 1:L 
and references therein), where a quantum paramagnetic 
phase near J2 ~ 0.5Ji is observed the nature of which 
is still under discussion. Another canonical model is 
the Shastry-Sutherland antiferromagnet introduced in 
the eighties^, which has special arrangement of frustrat- 
ing next-nearest-neighbor J2 bonds on the square lattice, 
cf. Fig^ We note that for bonds of equal strength, 
i.e., J\ — J2, the Shastry-Sutherland model is equiva- 
lent to a Heisenberg model on one of the eleven uniform 
Archimedean lattices^. Although the initial motivation 
to study this special frustrated square-lattice antiferro- 
magnet is related to the existence of a simple singlet- 
product eigenstate (which becomes the ground state (GS) 
for strong frustration), the renewed interest in the last 
years was stimulated by the discovering of the new quan- 
tum phase in SrCu(B0 3 )2^ which can be understood 
in terms of the Shastry-Sutherland model. Although the 
GS of this model in the limit of small frustration J2 and 



large J2 is well understood, the GS phase at moderate J2 
is still a matter of discussion. 

In this paper, we study the GS phase diagram for 
spin half Shastry-Sutherland model using a high-order 
coupled cluster treatment. The coupled cluster method 
(CCM) has previously been applied to various quantum 
spin systems with much succesa 8 ' 15 i 16 i 17 i 18 ' 19 i 20 i 21 i 22 i 23 . 
We mention that one particular advantage of this ap- 
proach consists in applicability to strongly frustrated 
quantum spin systems in any dimension, where some 
other methods, such as, e.g., the quantum Monte Carlo 
method fail. 



II. THE MODEL 

The Shastry-Sutherland model is a spin-i Heisenberg 
model on a square lattice with antiferromagnetic nearest- 
neighbor bonds Ji and with one antiferromagnetic diag- 
onal bond J2 in each second square (see Fig0. It is 
described by the Hamiltonian 

h = Ji Si ■ Sj + j 2 Si ■ (1) 

where the operators Sj represent spin-half operators, i.e., 
Sj = s(s + 1) with s = 1/2. The sums over and 
{i, k} run over all nearest-neighbor bonds and over some 
of the next-nearest-neighbor bonds according to the pat- 
tern shown in Fig. ^ Due to the special arrangement of 
the J2 bonds the unit cell contains four sites. Therefore 
it is convenient to split the square lattice into four equiv- 
alent sublattices A, C, B and D as shown in Fig.JTJ. It 
what follows we set J\ = 1 and consider J2 > as the 
parameter of the model. 
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FIG. 1: Illustration of the classical spiral state for Shastry- 
Sutherland model of Eq. with nearest-neighbor bonds 
Ji (solid lines) and next-nearest-neighbor bonds J2 (dashed 
lines). The spin orientations at A, C and B, D lattice sites 
are defined by the angles 9 — ncj) and 9 = ntfi + n, respectively, 
where n = 0, 1, 2, . . . , and <f> is the characteristic angle of the 
spiral state. The state is shown for <f> = tv/6 and n = 0, 1 ... 5. 



The classical (i.e., s — > 00) GS of the Shastry- 
Sutherland model is the collinear Neel state for J2/J1 < 
1, but a noncollinear spiral state for J 2 / J\ > 1 (see FigQ] 
and Rcfs. with a characteristic pitch angle <j) given 

by 







arccos(— J1/J2) 



■h < Ji 
J 2 > J 1 



(2) 



We note that for = the spiral state becomes the 
collinear Neel state classically. The transition from the 
collinear Neel to noncollinear spiral state is of second 
order and takes place at J2/J1 = 1- We note further that 
there are only two different angles between interacting 
spins, namely, 4> + ir for the J\ couplings and —2(f) for the 
J 2 couplings. 

The quantum s — 1/2 version of the model has been 
treated previously by various methods like Schwinger bo- 
son mean-field theory^, exact diagonalizatioi 



.14,26 



,25,27.28,29 



expansions*************, renormalization group and also by 
a gauge-theoretical approach 31 . A recent review can be 
found in Ref. |U From these studies one knows that for 
small J2 < J\ the physics of the quantum model is similar 
to that of the classical model, i.e., we have semi-classical 
Neel order. Furthermore, one knows already from the 
early work of Shastry and Sutherland^ that for large J2 
the quantum GS is a rotationally invariant product state 
of local pair singlets (so-called orthogonal-dimer state) 

\*) di mer = Y[ {h]] J U)\ ij) - I U)\ h)V ^ > where 

i and j correspond to those sites which cover the J2 
bonds. The energy per site of this orthogonal-dimer state 



is Edimer/N = — 3J2/8. It becomes the GS at around 
Jf w (1.44... 1.49) Ji (see Table 2 in Ref. HI). Note 
that such an orthogonal-dimer state can be observed also 
in corresponding one-dimensional and three-dimensional 
m0( j e i s 33 4 34 i 35 rpj^g na ^ ure f transition between the 

semi-classical Neel state and the orthogonal-dimer phase 
is still a matter of controversial discussion. In the region 
1.2Ji < J 2 < 1.45Ji the main question is whether the 
system has an intermediate phase. A direct transition be- 
tween the Neel phase and the orthogonal-di mer phase is 
favored in Refs. Ill27ll2l8(i whereas in Refs. 1241251261311 
the existence of an intermediate phase is found. How- 
ever, concerning the nature of this intermediate phase 
controversial results are reported, as candidates for the 
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or a 



intermediate phase are quantum spiral phase, 
plaquette or columnar singlet phase s 25 ! 26 discussed. 

To contribute to the solution of this open problem the 
CCM is an appropriate method, since it is one of the 
methods which can deal with spiral phases in quantum 
spin modelsi 19 i 21 i 36 



III. THE COUPLED CLUSTER METHOD 

The CCM formalism is now briefly considered, al- 
th ough t he interested reader is referred to Refs. ITalTsl 
I19H2O1I22I for further details. The starting point for the 
CCM calculation is the choice of a normalized reference 
or model state |<E>). For spin systems, an appropriate 
choice for the CCM model state |$) is often a classical 
spin state, in which the most general situation is that 
each spin can point in an arbitrary direction. In order 
to treat the Shastry-Sutherland model using the CCM, 
we choose the Neel state and the spiral state in Fig.Q 
to be our model states. We note that we do not choose 
the classical result for the pitch angle <j> but we consider 
it rather as a free parameter in the CCM calculation. 

To treat each site equivalently we perform a rotation 
of the local axis of the spins such that all spins in the 
reference state align in the same direction, namely along 
the negative z axis, such that we have | < I ) ) = ||) ||) \[) ■ ■ ■ ■ 
We define a set of multi-spin creation operators Cf = 

s t ) S t S l ' S -r S l S m j • • • • 

The choice of the Cf ensures that ($|C/ = = C/|$), 
where Cj is the Hermitian adjoint of CJ . 

In order to make the spin Si to be aligned along the 
negative z axis one has to perform a rotation of the re- 
spective spin by an appropriate angle <5j. This rotation 
is equivalent to the canonical transformations, 

sf — cos <5jsf -I- sin 5isf 

4 = 4 

s| = — sin 5isf + cos ^sf. (3) 
Using this transformation the Hamiltonian is then 
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rewritten as 

N 

H = 3 X ^ (2 sin( Pi,j [ § t § j - § i s t + S i §Z i ~ S i S 7^ 

(id) 

+ CQB<p i , j sZs* + j(cOSy>i,j + l)[sfsj + §-§ + } 

+i(cosy>ij - l)[s 4 + s+ + s^sj]) 

N 

+ cos^ a s?s£ + -{coscp itk + l)[sfs^ +s~s+] 

+ i(cos^ fc - l)[sfs+ +s-S-]y (4) 

where the angles ifiij = Sj — Si, ip^k = 5k — Si between 
two nearest-neighbor and next-nearest-neighbor spins are 
ifiij =w + (j), tfi.k — —20, respectively, and s ± = s x ±is v 
are spin raising and spin lowering operators. 

The ket and bra GS's |*) and (*| of H are 
parametrised within the CCM as follows: 

H\$) = E\V) ; (W\H = E{$\ ; 
|*>=e s |$) ; S = Y,SiC+ ; 

(*| = ($|5e- s ; S=l + J2SiCT- (5) 

IjiO 

The correlation operators S and S contain the corre- 
lation coefficients Si and Si which have to be deter- 
mined. Using the Schrodinger equation, H\^>) = E\^), 
we can now write the GS energy as E — (Q\e~ s He s \&). 
After the notational rotation of the local axes of the 
quantum spins, the sublattice magnetization is given by 
M = -l/tf£f<*l*?l*>. 

To find the ket-state and bra-state correlation coeffi- 
cients Si and Si we require that the expectation value 
H = i^\H\^) is a minimum with respect to Si and Si, 
such that the CCM ket-state and bra-state equations are 
given by: 

(<S>\CYe- s He s \<S>) =0 V/ ^ 
($\Se- s [H, C+]e s |$> =0 V/ ^ 0. 

(6) 

The CCM formalism is exact if we take into account all 
possible multispin configurations in the correlation oper- 
ators S and S, which is, however, in general impossible 
for a quantum many-body model. Hence, it is neces- 
sary to use approximation schemes in order to truncate 
the expansion of S and S in the Eqs. JHJ in any practi- 
cal calculation. The most common scheme is the LSUBn 
scheme, where we include only n or fewer correlated spins 
in all configurations (or lattice animals in the language 
of graph theory) which span a range of no more than n 
adjacent (contiguous) lattice sites. 



TABLE I: Number of fundamental GS configurations of the 
LSUBn approximation for the Shastry-Sutherland model us- 
ing the Neel state (<j> — 0) and the spiral state (<f> 7^ 0) as the 
CCM reference state. 



LSUBn 


Neel state: <f> = 


spiral state: cj> 7^ 


2 


1 


12 


4 


35 


248 


6 


794 


6184 


8 


20892 


166212 



To find all possible fundamental configurations which 
are different under the point and space group symmetries 
of both the lattice and the Hamiltonian, we use the lattice 
symmetries. The numbers of fundamental configurations 
may be further reduced by the use of additional conser- 
vation laws. For example, in the case of the Neel state 
(4> = 0), the Hamiltonian of Eq. commutes with the 
total uniform magnetization, = ^2 k s| (the sum on k 
runs over all lattice sites), the GS lies in the sf, = sub- 
space, and hence we exclude configuration with an odd 
number of spins or with unequal numbers of spins on the 
two equivalent sublattices. For the spiral state we cannot 
apply this property because it is not an eigenstate of Sj,. 
We calculate the fundamental configurations numerically, 
and the results of the numbers of LSUBn configurations 
for n < 8 are given in Table I. By using parallel com- 
puting we are able to solve the 20892 equations of the 
CCM-LSUB8 approximation for the Neel reference state. 
However, for the spiral state the current limitations of 
computer power allow then solution of the CCM equa- 
tions up to LSUB6, only. 

Since the LSUBn approximation becomes exact in the 
limit n — > 00, it is useful to extrapolate the 'raw' LSUB?i 
results to the limit n — > 00. Although an exact scal- 
ing theory for the LSUBn results is not known, there is 
some empirical experiencei&i&Si how the physical quan- 
tities for antiferromagnetic spin models scale with n. As 
stated above for the Neel reference state we are able to 
calculate the GS energy E and the sublattice magnetiza- 
tion M within LSUBn up to n = 8. In order to obtain 
more accurate results for the GS energy, we now employ 
a scaling Iawi2i2£ in order to extrapolate our results in 
the limit m — > 00, where 

E(n) = a + ai\ + a 2 ( \ J . (7) 

We use CCM results for n = 4, 6, 8 in order to carry 
out these extrapolations. 20 We find, however, that other 
scaling laws proposed in the literature yield very simi- 
lar results for the energy. In the Neel ordered phase we 
utilise— a scaling law with leading power 1/n, i.e., 

M(n) = b Q + b 1 -+b 2 (-) . (8) 
n \n J 

We find that this prescription again leads to reasonable 
results 20 . However, applying this scaling rule to systems 
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showing an order-disorder transition at zero temperature 
this kind of scaling tends to overestimate the magnetic 
order and yields too large critical values for the exchange 
parameter driving the transitioni2*2£. The reason for that 
might consist in the change of the scaling near a critical 
point. Hence in addition to the scaling rule (JHJ) we also 
use a leading 'power-law' scaling^, given by 

M(n) = c + cJ^\ . (9) 

The leading exponent ci is determined directly from the 
LSUBn data. 



IV. RESULTS 
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We start with the discussion of the onset of the spiral 
phase in the quantum model. We calculate the GS en- 
ergy as a function of Ji using as reference state a spiral 
state as sketched in Fig. ^ As quantum fluctuations 
may lead to a "quantum" pitch angle that is different 
from the the classical case, we consider the pitch angle 
in the reference state as a free parameter. We then de- 
termine the "quantum" pitch angle <j) qu by minimizing 
£'LSUBm(0) with respect to </> in each order n. As for 



FIG. 3: The "quantum" pitch angle 



as a function of 



E/N 




-0.5285 
-0.529 



FIG. 2: Ground-state energy versus the pitch angle <j> within 
CCM-LSUB4 approximation for different values of J2 in the 
range 1.55 < J 2 < 1.59. 

the classical model for small J2 the energy i^LSUBm (</>) 
has its minimum at <j) qu = 0, i.e., the quantum GS is the 
semi-classical collinear Neel state. Contrary to the classi- 
cal case, this collinear quantum state can survive into the 
region Ji > J%, where classically it is already unstable. 
This effect is known as order from disordei^^ and is 
widely observed in quantum spin systems, see, e.g., Refs. 
Il9l3a . For frustrating couplings J2 ^ I.5J2 apart from 
the minimum at 4> = a second minimum at a finite 
<p > emerges, which becomes the global minimum for 
strong enough J2. This scenario illustrated in Fig. [3 is 



J2 calculated within CCM-LSUBn approximation with n = 
2,4,6. 



typical for a first-order transition, i.e., we find indications 
that quantum fluctuations may change the nature of the 
phase transition between the the collinear Neel phase to 
the noncollincar spiral phase from a second-order classi- 
cal transition to a first-order quantum transition. Note 
that a similar situation can be found in other frustrated 



spin system: 



,19.21 



The "quantum" pitch angle (j) qu , where 



EhSVBm(<P) has its global minimum, is shown in Fig. [3] 
(j> qu shows a typical jump from <p qu = to a finite value. 
Our data clearly indicate that the quantum noncollinear 
spiral phase has lower energy than the collinear phase 
only for strong frustration J2 > 1.5Ji. 
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FIG. 4: The energy of (i) the collinear quantum ground state 
as function of J2 obtained by CCM-LSUBn with n = 4, 6, 8 
and its extrapolated value to n — > 00, see Eq. 0, and (ii) of 
the orthogonal-dimer state. 



Next we compare the energy of the orthogonal-dimer 
state \ty)dimer and the energy of the collinear quantum 
ground state (i.e. the reference state |<&o) is the Neel 
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state), see Fig 0] We can postpone the discussion of 
the question whether that quantum ground state pos- 
sesses Neel LRO or not, since it is possible (starting 
from the Neel reference state) to calculate the energy 
up to high accuracy even in a parameter regime where 
the Neel order breaks down due to quantum fluctua- 
tions, i.e. for a magnetically disordered state, see e.g. 
Refs. I8EHI9II21II2I2I Our results demonstrate, that 
the orthogonal-dimer state has lower energy than the 
collinear state for J 2 > 1.477 J%. \^)dimer remains the 
state of lowest energy also in the region where the non- 
collinear spiral state has lower energy than the collinear 
phase. We conclude that there is no intermediate spi- 
ral phase in the quantum model. Our estimate of the 
critical value Jf* = 1.477 Ji where the transition to the 
orthogonal-dimer phase takes place is in g ood agreement 
with other results, cf. Table 2 in Ref. l32t 

So far we have discussed mainly the energy of com- 
peting GS phases. The last question we would like to 
discuss is the question of the stability of the Neel LRO 
in the frustrated regime. For that we calculate the or- 
der parameter (sublattice magnetization) M within the 
LSUBn approximation scheme up to n = 8 and extrapo- 
late to n — > 00 using two variants of extrapolation as de- 
scribed in Sect. IIIII The results are shown in Fig. O The 
extrapolated data clearly demonstrate that the LRO van- 
ishes before the orthogonal-dimer state becomes the GS. 
The transition from Neel LRO to magnetic disorder is of 
second order. Hence we come to the second important 
statement that there exists an intermediate magnetically 
disordered phase. Within the used CCM scheme starting 
from the Neel reference state we are not able to discuss 
the nature of the magnetically disordered state preceding 
the orthogonal-dimer state. Though there are some first 
attempts to develop a CCM formalism for magnetically 
disordered valence bond phases^, a high level of approx- 
imation is reached currently only starting with Neel or 
spiral reference states. 

Obviously, the critical value where Jf the Neel LRO 
breaks down depends on the used extrapolation formula. 
The extrapolation according to Eq. JSJl leads accurate re- 
sults for M in the unfrustrated (J2 = 0) square-lattice 
limit and yields Jf ~ 1.39Ji. As discussed in Sect. IIIII 
this extrapolation scheme tends to overestimate the re- 
gion of magnetic LRO and indeed the value Jf/Ji = 1-39 
is significantly larger than the corresponding value cal- 
culated by series expansion, see Table 2 in Ref. The 
extrapolation according to Eq. (JSJ with a variable ex- 



ponent C2 is less accurate in the unfrustrated limit but 
it seems ro be more appropriate to find the position of 
the critical point Jf , since the scaling behavior might be 
changed at the critical point. We get Jf ~ 1.14Ji which 
fits well to the corresponding value calculated by series 
expansion. 
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FIG. 5: Sublattice magnetization M versus J2 obtained by 
CCM-LSUBn with n — 4, 6, 8 and its extrapolated values to 
n — > 00 using two different extrapolation schemes, namely 
according to Eq. JH) (extrapol 1) and to Eq. © (extrapol 2). 



CONCLUSIONS 



We have studied the GS phase diagram of the spin 
half Shastry-Sutherland amtiferromagnet making use of 
high-order coupled cluster calculations. Comparing the 
energies of competing Neel , spiral and orthogonal-dimer 
phases we can rule out the existence of a noncollinear 
spiral phase. Considering the Neel order parameter we 
find that the semi-classical Neel long-range order disap- 
pears before the orthogonal-dimer phase sets in. Hence 
we conclude that the Neel phase and the dimer phase 
are separated by a magnetically disordered intermediate 
phase. 
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